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The paper deals with the Weyl equation which is the massless Dirac equation. We study the 
Weyl equation in the stationary setting, i.e. when the spinor field oscillates harmonically in time. 
We suggest a new geometric interpretation of the stationary Weyl equation, one which does not 
require the use of spinors, Pauli matrices or covariant differentiation. We think of our 3-dimensional 
space as an elastic continuum and assume that material points of this continuum can experience no 
displacements, only rotations. This framework is a special case of the Cosserat theory of elasticity. 
Rotations of material points of the space continuum are described mathematically by attaching 
to each geometric point an orthonormal basis which gives a field of orthonormal bases called the 
coframe. As the dynamical variables (unknowns) of our theory we choose the coframe and a density. 
We choose a particular potential energy which is conformally invariant and then incorporate time 
into our action in the standard Newtonian way, by subtracting kinetic energy. The main result of 
our paper is the theorem stating that in the stationary setting our model is equivalent to a pair of 
Weyl equations. The crucial element of the proof is the observation that our Lagrangian admits a 
factorization. 
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I. MAIN RESULT 

Throughout this paper we work on a 3-manifold M 
equipped with local coordinates x", a — 1, 2, 3, and pre- 
scribed positive metric gap which does not depend on 
time. We extend the Riemannian 3-manifold M to a 
Lorentzian 4-manifold M x M by adding the time coor- 
dinate G M. The metric on R x M is defined as 



equation is 
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Here and further on we use bold type for extended quan- 
tities. Say, the use of bold type in tensor indices a,f3 
appearing in the LHS of formula ([1]) indicates that these 
run through the values 0,1,2,3, whereas the use of nor- 
mal type in tensor indices a, /? appearing in the RHS of 
formula ([ij indicates that these run through the values 
1,2,3. 

All constructions presented in the paper are local so 
we do not make a priori assumptions on the geometric 
structure of {M,g}. 

The Weyl equation is the massless Dirac equation. It 
is the accepted mathematical model for a massless neu- 
trino field. The dynamical variable (unknown quantity) 
in the Weyl equation is a two-component complex- valued 
spinor field ^ which is a function of time g R and lo- 
cal coordinates x" on M. The explicit form of the Weyl 
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Here the a are Pauli matrices, da is the time deriva- 
tive and Vet is the covariant spatial derivative, see Ap- 
pendix|B]for details. Summation in ([2]) is carried out over 
the tensor index a = 1,2,3 as well as over the spinor in- 
dex b = 1,2. The use of the partial derivative do = d/dx'^ 
in equation ([2|) is justified by the fact that the time co- 
ordinate x^ is fixed and we allow only changes of coordi- 



nates {x^ , 



^) which do not depend on a;°. 



We see that the Weyl equation ([2]) is a system of two 
(d = 1,2) complex linear partial differential equations 
on the 4-manifold R x Af for two complex unknowns 
6 = 1,2. The two choices of sign in ^ give two versions 
of the Weyl equation which differ by time reversal. Thus, 
we have a pair of Weyl equations. 

We will be interested in spinor fields of the form 



where 



^(x^a;^a;^x3) = e-'P'>'=%^{x\x^ ,x^) (3) 



Po ^ (4) 



is a real number. Substituting ^ into ^ we get the 
equation 







(5) 



which we shall call the stationary Weyl equation. The 
difference between equations ([2]) and (O is that in equa- 
tion ((2]) the spinor field ^ "lives" on the Lorentzian 4- 
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1] "lives" on the Riemannian 3-manifold M. Thus, the 
stationary Weyl equation is the Weyl equation with time 
separated out. The relationship between equations ^ 
and ^ is the same as between the wave equation and 
the Helmholtz equation. 

The stationary Weyl equation ^ is the object of study 
of this paper. We separated out time to simplify the prob- 
lem while retaining most of its essential features. Note 
also that this separation of variables has a clear physical 
meaning: the real number pq appearing in ^ and ([5]) is 
quantum mechanical energy. 

The aim of our paper is to show that the stationary 
Weyl equation ([5]) can be reformulated in an alterna- 
tive (but mathematically equivalent) way using instead 
of a spinor field a different set of dynamical variables. 
Namely, we view our 3-manifold M as an elastic contin- 
uum whose material points can experience no displace- 
ments, only rotations, with rotations of different mate- 
rial points being totally independent. The idea of ro- 
tating material points may seem exotic, however it has 
long been accepted in continuum mechanics within the 
Cosserat theory of elasticity [ij. This idea also lies at 
the heart of the theory of teleparallelism (= absolute 
parallelism = fernparallelismus), a subject promoted by 
A. Einstein and E. Cartan 0-13 . See Section E] for more 
details. 

Rotations of material points of the 3-dimensional elas- 
tic continuum are described mathematically by attach- 
ing to each geometric point of the manifold M an or- 
thonormal basis, which gives a field of orthonormal bases 
called the frame or coframe, depending on whether one 
prefers dealing with vectors or covectors. Our mathemat- 
ical model will be built on the basis of exterior calculus 
(no need for covariant derivatives) so for us it will be 
more natural to use the coframe. 

The coframe is a triple of orthonormal covector fields 
'd^ , j = 1, 2, 3, on the 3-manifold M. Each covector field 
'd^ can be written more explicitly as i9^q, where the ten- 
sor index a — 1,2,3 enumerates the components. The 
orthonormality condition for the coframe can be repre- 
sented as a single tensor identity 

g = ® (6) 

where Sjk is the Kronecker delta. For the sake of clarity 
we repeat formula ([6]) giving tensor indices explicitly and 
performing summation over frame indices explicitly: 

where a and /3 run through the values 1,2,3. We view 
the identity (jH) as a kinematic constraint: the metric 
g is given (prescribed) and the coframe elements -d^ are 
chosen so that they satisfy ^ , which leaves us with three 
real degrees of freedom at every point of M. 

As dynamical variables in our model we choose the 
coframe i? and a positive density p. Our coframe and 
density are functions of local coordinates {x^,x'^,x^) on 
M as well as of time a;°. 



At a physical level, making the density p a dynamical 
variable means that we view our continuum more like a 
fiuid rather than a solid. In other words, we allow the 
material to redistribute itself so that it finds its equilib- 
rium distribution. 

Note that the total number of real dynamical degrees of 
freedom contained in the coframe and positive density 
p is four, exactly as in a two-component complex-valued 
spinor field Moreover, it is known (see Appendix [C]) 
that a coframe •& and a (positive) density p are geometri- 
cally equivalent to a nonvanishing spinor field ^ modulo 
the sign of 

The crucial element in our construction is the choice of 
potential energy. It is known that in the purely rotational 
setting the potential energy of a physically linear elas- 
tic continuum contains three quadratic terms, with three 
real parameters (elastic moduli) as factors. The three 
quadratic terms in potential energy correspond to the 
three irreducible pieces of torsion. It is not a priori clear 
what the elastic moduli of "world aether" are. We choose 
a potential energy which feels only one piece of torsion, 
axial, and is, moreover, conformally invariant, i.e. which 
is invariant under the rescaling of the 3-dimensional met- 
ric g by an arbitrary positive scalar function. This leaves 
us with a unique, up to rescaling by a positive constant, 
formula ([T^ for potential energy. 

After the potential energy is chosen the remainder of 
our construction is straightforward. We incorporate time 
into our action in the standard Newtonian way, by sub- 
tracting kinetic energy. This gives us the Lagrangian den- 
sity As we are interested in comparing our math- 
ematical model with the Weyl equation, we perform a 
change of dynamical variables and switch from coframe 
i9 and density p to a spinor field ^. Our Lagrangian den- 
sity now takes the form ([^0)1 . We write down the field 
equation (Euler-Lagrange equation) for our Lagrangian 
density and observe that time separates out if we seek sta- 
tionary solutions Q ; this separation of variables is highly 
nontrivial because our field equation is nonlinear. After 
separation of variables our Lagrangian density takes the 
stationary form ([25)1 . 

The main result of our paper is the following 

Theorem 1 A nonvanishing time-independent spinor 
field rj is a solution of the field equation for our station- 
ary Lagrangian density ^23^ if and only if it is a solution 
of one of the two stationary Weyl equations 

Theorem [T] provides an elementary, in terms of New- 
tonian mechanics and elasticity theory, interpretation of 
the stationary Weyl equation. This interpretation is ge- 
ometrically much simpler than the traditional one as the 
mathematical description of our model does not require 
the use of spinors, Pauli matrices or covariant differenti- 
ation. 

The only technical assumption contained in the state- 
ment of Theorem [Hand its proof is that the density does 
not vanish which is equivalent to the spinor field not van- 
ishing. At the moment we do not know how to drop this 
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technical assumption. We can only remark that generi- 
cally one would not expect a spinor field 77 "living" on a 
3-manifold to vanish as this would mean satisfying four 
real equations Re 77^ —Imi]^ = Re 77^ = Imrj^ ~ having 
at our disposal only three real variables , a = 1,2,3. 

The crucial element of the proof of Theorem [T] is the 
observation that our Lagrangian density admits factor- 
ization, see formula (1311) . Thus, our argument is similar 
to the original argument of Dirac, the difference being 
that we factorize the Lagrangian whereas Dirac factor- 
ized the field equation (Klein-Gordon equation). In our 
model factorizing the field equation is impossible because 
the equation is nonlinear. 

Our paper has the following structure. In Section |ll] 
we describe our mathematical model, in Section IIIII we 
switch to the language of spinors, in Section ITVl we sep- 
arate out time and in Section |V] we factorize our La- 
grangian. Section I VII contains the proof of Theorem [T] 
In Section IVIII we analyze plane wave solutions and in 
Section IVIIII we give a relativistic representation of our 
Lagrangian. The concluding discussion is presented in 
Section nil 

Our paper also has five appendices dealing with is- 
sues of a more technical nature: Appendix |X] describes 
our general notation. Appendix |B] describes our spinor 
notation. Appendix [C] gives the correspondence between 
coframes and spinors. Appendix |D] gives the spinor rep- 
resentation of axial torsion and angular velocity and Ap- 
pendix [E] deals with a toy model which illustrates how 
a particular class of nonlinear second order differential 
equations reduces to pairs of linear first order equations. 

II. OUR MODEL 

In this section we describe in detail our mathematical 
model. At a basic level it was already sketched out in 
Section [D 

We need to write down the potential energy of a de- 
formed Cosserat continuum. The natural measure of de- 
formations caused by rotations of material points is the 
torsion tensor defined by the explicit formula 

T Syk-ff^ ® dd'' (7) 

where d denotes the exterior derivative. Here "tor- 
sion" means "torsion of the teleparallel connection" with 
"teleparallel connection" defined by the condition that 
the covariant derivative of each coframe element is 
zero; see Appendix A of Q for a concise exposition. 

Our construction of potential energy follows the logic 
of classical linear elasticity 0, the only difference be- 
ing that instead of a rank 2 tensor (strain) we deal with 
a rank 3 tensor (torsion). The logic of classical linear 
elasticity dictates that we must first decompose our mea- 
sure of deformation (torsion) into irreducible pieces, with 
irreducibility understood in terms of invariance under 
changes of local coordinates preserving the metric gap 
at a given point P S Af. It is known [7| that torsion has 



three irreducible pieces labeled by the adjectives axial, 
vector and tensor. (Vector torsion is sometimes called 
trace torsion.) The general formula for the potential 
energy of a homogeneous isotropic linear elastic mate- 
rial contains squares of all irreducible pieces with some 
constant coefficients in front. Thus, the general formula 
for potential energy should contain three free parameters 
(elastic moduli). 

We, however, choose to construct our potential energy 
using only one piece of torsion, namely, the axial piece 
given by the explicit formula 

T^'^ i,5jfei?J' Adi?^ (8) 

Comparing ^ with ([T]) we see that axial torsion has a 
very simple meaning: it is the totally antisymmetric part 
of the torsion tensor (T is antisymmetric only in the last 
pair of indices whereas T^^ is antisymmetric in all three) . 
In other words, T^^ is a 3-form. 

We chose the axial piece of torsion because it has two 
remarkable properties. 

• The definition of axial torsion ([8|) is very simple in 
that it does not involve the metric. In a sense, axial 
torsion (3-form) is an analogue of the electromag- 
netic field tensor (2-form) from Maxwell's theory. 

• Axial torsion possesses the property of conformal 
covariance, i.e. scales nicely under conformal rescal- 
ings of the metric. Indeed, it is easy to see that if 
we rescale our coframe as 

^ e^§'^ (9) 

where /i : M — >■ M is an arbitrary scalar function, 
then our metric scales as 

5a/3 ^ e^^g^p (10) 
and axial torsion scales as 

without the derivatives of h appearing. The fact 
that axial torsion is conformally covariant was 
previously observed by Yu. N. Obukhov [1] and 
J. M. Nester @. 

We take the potential energy of our continuum to be 
P(a;") := / WT^^'fpdxHx^dx^. (12) 

J M 

It is easy to see that the potential energy is con- 
formally invariant: it does not change if we rescale our 
coframe as ^ and our density as 

p ^ e^'^p. (13) 

This follows from formulas dTTI) . ([TU)) and 

||rTiax||2 ^ rpax T-iax ^OiK BX 'yu 

II — -^^^aPi^ ii\ti9 9 9 ■ 
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We take the kinetic energy of our continuum to be 

K{x°) := / W^fpdx^dx^dx^ (14) 

where i? is the 2-form 

^ y^k^' A do^" (15) 

(compare with ([5])). The 2-form (IT5|) can, of course, be 
written as 

2 

i?=-*a; (16) 

where 

^ := ^ * A aoz?") (17) 

is the (pseudo)vector of angular velocity. Hence, 
is the standard expression for the kinetic energy of a 
homogeneous isotropic Cosserat continuum. In writing 
formula (|14p we assumed homogeneity (properties of the 
material are the same at all points of the manifold M) 
and isotropy (properties of the material are invariant un- 
der rotations of the local coordinate system) . We think 
of each material point as a uniform ball possessing a mo- 
ment of inertia and without a preferred axis of rotation. 

We now combine the potential energy ([T2|) and kinetic 
energy (fH)) to form the action (variational functional) of 
our dynamic problem: 

S{d,p) J (P(a;0) - if (a;°)) dx° 

= [ L{i},p)dx°dx^dx^dx^ (18) 

JrxM 

where 

L(^,p):=(||r-|p-||^||> (19) 

is our Lagrangian density. Note that our construction of 
the action out of potential and kinetic energies is 
Newtonian (compare with classical elasticity or even the 
harmonic oscillator in classical mechanics). 

Our field equations (Euler-Lagrange equations) are ob- 
tained by varying the action with respect to the 
coframe and density p. Varying with respect to the den- 
sity p is easy: this gives the field equation HT'^'^lp ~ 
which is equivalent to L(j!),p) = 0. Varying with re- 
spect to the coframe is more difficult because we have 
to maintain the kinematic constraint recall that the 
metric is assumed to be prescribed (fixed). 

A technique for varying the coframe with kinematic 
constraint ^ was described in Appendix B of Q. We, 
however, do not write down the field equations for the 
Lagrangian density (|19l) explicitly. We note only that 
they are highly nonlinear and do not appear to bear any 
resemblance to the linear Weyl equation ([2]). 



Remark 1 The 3-form T'^^ and 2-form are invari- 
ant under rigid rotations of the coframe, i. e. under spe- 
cial orthogonal transformations HAS]) with constant k- 
Hence, our Lagrangian density U9\) is invariant under 
rigid rotations of the coframe and, accordingly, solutions 
of our field equations whose coframes differ by a rigid 
rotation can be collected into equivalence classes. Fur- 
ther on we view coframes differing by a rigid rotation as 
equivalent. 

III. SWITCHING TO THE LANGUAGE OF 
SPINORS 

As pointed out in the previous section, varying the 
coframe subject to the kinematic constraint © is not an 
easy task. This technical difficulty can be overcome by 
switching to a different dynamical variable. Namely, it is 
known, see Appendix [Cl that in dimension 3 a coframe 
"i? and a (positive) density p are equivalent to a nonva- 
nishing spinor field ^ modulo the sign of ^. The great 
advantage of switching to a spinor field ^ is that there 
are no kinematic constraints on its components, so the 
derivation of field equations becomes absolutely straight- 
forward. 

We now need to substitute formulas (jCl[) . (jC3p and 
([Cl]) into ^ and UHl) to get explicit expressions for T"^^ 
and 1? in terms of the spinor field The results are pre- 
sented in Appendix |DJ Namely, formula (jPlI) gives the 
spinor representation of the 3-form T^^ whereas formu- 
las (|D2p and give the spinor representation of the 
2-form {). We also know the spinor representation for 
our density p, see formulas (jCip and (|C2p . Substitut- 
ing all these into formula (IT5|) we arrive at the following 
self-contained explicit spinor representation of our La- 
grangian density 

- Ue<yaabdoe - e'^adtaof)|p) Vd^. (20) 

Here and further on we write our Lagrangian density 
and our action as L(^) and S{£,) rather than L{d,p) and 
S{'d, p), thus indicating that we have switched to spinors. 
The nonvanishing spinor field ^ is the new dynamical 
variable and it will be varied without any constraints. 

Straightforward calculations show that the field equa- 
tion for our Lagrangian density (j20p is 

Ai 

- -[{*T-^)cj^abVo.e + <T'^dbVc,{{*T-^)e)) 

y 

~ p-^LaoabS!" = Q (21) 

where the geometric quantities *T'*'^, oj, p and L are ex- 
pressed via the spinor field ^ in accordance with formu- 
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las dm]), dm]), dCT]), dC2]) and ^Q. The LHS of equa- 
tion is the spinor field appearing in the formula 
for the variation of the action (flSl): 



JrxM 

We shall refer to equation dHJ as the dynamic field 
equation, with "dynamic" indicating that it contains the 
time derivative do- 



IV. SEPARATING OUT TIME 



Our dynamic field equation dS]) is highly nonlinear 
and one does expect it to admit separation of variables. 
Nevertheless, we seek solutions of the form dS])- Substi- 
tuting formula d3 into formulas dE]), dHl]), dCl]), dC2|) 
and (|20|) and using the identity dBSP we get 



*J = — , [22) 



2poT]''aadbV 



d 



(23) 



(24) 
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9r]^(JocdV 



[pot^^dbrtf^^A^g. (25) 



Note that the geometric quantities d^^ ~ d^5|) do not de- 
pend on time x^ , which simplifies the next step: substi- 
tuting dSJ into our dynamic field equation (|2ip , using the 
identity (IBSp and dividing through by the common factor 

4? 

- y ((*T-")'T"d6V„77'' + a«a6Va((*T'^^)77'')) 
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-p-iLaodbr;'' = 0. (26) 



The remarkable feature of formulas dHJ-dSH]) is that 
they do not contain dependence on time . Thus, we 
have shown that our dynamic field equation (I2ip admits 
separation of variables, i.e. one can seek solutions in the 
form dS]). 

We shall refer to equation d^ as the stationary field 
equation, with "stationary" indicating that time x^ has 
been separated out. 

Consider now the action 

S{7j):= [ L{7j)dx^dx^dx^ (27) 



where L^ij) is our "stationary" Lagrangian density d25l) . 
It is easy to see that our stationary field equation dM]) 
is the Euler-Lagrange equation for our "stationary" ac- 
tion (l?f)) . 

In the remainder of the paper we do not use the ex- 
plicit form of the stationary field equation I^E^, dealing 
only with the stationary Lagrangian density (I25p and the 
stationary action d^7)) . We needed the explicit form of 
field equations, dynamic and stationary, only to justify 
separation of variables. 

It appears that the underlying group-theoretic reason 
for our nonlinear dynamic field equation d2ip admitting 
separation of variables is the fact that our model is U(l)- 
invariant, i.e. it is invariant under the multiplication of 
the spinor field ^ by a complex constant of modulus 1. 
Hence, it is feasible that one could have performed the 
separation of variables argument without even writing 
down the explicit form of field equations. 

We give for reference a more compact representation of 
our stationary Lagrangian density dUI in terms of axial 
torsion T'^^ (see formula d22)) ) and density p (see for- 
mula dM])): 



L(,?)=(||r-|p_^p5jp 



(28) 



Of course, formula dH]) is our original formula dlU) 
with time separated out. The choice of dynamical vari- 
ables in the stationary Lagrangian density l\28\i is up 
to the user: one can either use the time-independent 
spinor field ry or, equivalently, the corresponding time- 
independent coframe and time-independent density (the 
latter are related to rj by formulas dClJ-dC4| with ^ re- 
placed by rj). The important thing is that now our dy- 
namical variables are time-independent because we have 
separated out time. 

The fact that we use the same notation L both for the 
dynamic and stationary Lagrangian densities should not 
cause problems as in all subsequent sections, apart form 
Section IVIIIi we deal with the stationary case only. 



V. FACTORIZATION OF OUR LAGRANGIAN 



Put 



'-M^'^dbVo.ri' -Ti'a'^dbVo.t) 



Vdet g . (29) 



This is the Lagrangian density for the stationary Weyl 
equation dSJ. Formula d29p can be written in more com- 
pact form as 



(30) 



where ^T^^ is the Hodge dual of axial torsion, see for- 
mula (1^^ . and p is the density, see formula (1^^ . Com- 
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paring formulas (1281) and ((30|) we get 

32po L+(?7)L_(r;) 



9 L+(,7)-L_(r;) 



(31) 



Let us emphasize once again that throughout this pa- 
per we assume that the density p does not vanish, which 
is, of course, equivalent to the spinor field not vanishing. 
In view of formulas (pO)) and ^ in the stationary case 
the assumption p can be equivalently rewritten as 



L+(r;)^L_(r;) 



(32) 



so the denominator in is nonzero. 

Formula (PT|) is the centerpiece of our paper: it es- 
tablishes the connection between Cosserat elasticity and 
the Weyl equation. Moreover, the fact that the RHS of 
formula (j3ip contains a product of two Weyl Lagrangian 
densities shows that we are essentially following Dirac's 
factorization construction, the difference being that in 
the nonlinear setting we cannot factorize equations and 
have to settle for the next best thing — factorizing the 
Lagrangian. 



VI. PROOF OF THEOREM [T] 

Observe that the Lagrangian densities L± defined by 
formula (j29[) possess the property of scaling covariance: 



i±(e"r/) = e^'^L±(ry) 



(33) 



where h : M W is an arbitrary scalar function. In 
fact, the Lagrangian density of any formally selfadjoint 
(symmetric) linear first order partial differential operator 
has the scaling covariance property ([55)) . 

We claim that the statement of Theorem[l]follows from 
formulas (|3T|) and ([33| . The proof presented below is an 
abstract one and does not depend on the physical nature 
of the dynamical variable rj, the only requirement being 
that it is an element of a vector space so that scaling 
makes sense. 

Note that formulas (I5l]) and ([55)) imply that the La- 
grangian density L possesses the property of scaling co- 
variance, so all three of our Lagrangian densities, L, 
and L_, have this property. Note also that if 77 is a 
solution of the field equation for some Lagrangian den- 
sity C possessing the property of scaling covariance then 
£(ry) = 0. Indeed, let us perform a scaling variation of 
our dynamical variable 



rj i-^ + hr] 



(34) 



where ft, : M — > R is an arbitrary "small" scalar function 
with compact support. Then = SJ C{rj) = 2 J hC{rj) 
which holds for arbitrary h only if C{irj) = 0. 

In the remainder of the proof the variations of 77 are 
arbitrary and not necessarily of the scaling type ([34| . 

Suppose that 77 is a solution of the field equation for 
the Lagrangian density i+. [The case when 77 is a solu- 
tion of the field equation for the Lagrangian density L_ 



is handled similarly.] Then L+{r]) = and, in view of 
formula ([5^ . L^{'q) ^ 0. Varying -q we get 



L-iv) 



32po 

9 \J L+iTj) - L^ir,) 



L-iv) 



L+irj) - L^irj) 



32po 



32po 



so 



sim = ^siL+iv) 



(35) 



We assumed that r/ is a solution of the field equation 
for the Lagrangian density i+ so <5/L+ (77) = and for- 
mula pSI) implies that S J L[r}) = 0. As the latter is true 
for an arbitrary variation of 77 this means that 77 is a so- 
lution of the field equation for the Lagrangian density L. 

Suppose that 7/ is a solution of the field equation for the 
Lagrangian density L. Then L(r/) = and formula ([5T)) 
implies that either L^{rj) = or L^{rf) — 0; note that in 
view of ()32p we cannot have simultaneously L^{ri) = 
and L^{r]) — 0. Assume for definiteness that L+{r]) = 0. 
[The case when L-irf) = is handled similarly.] Vary- 
ing 77 and repeating the argument from the previous para- 
graph we arrive at ([55)) . We assumed that 77 is a solu- 
tion of the field equation for the Lagrangian density L so 
5 jL{ri) = and formula ([55)) implies that 5 J L+{r]) — 0. 
As the latter is true for an arbitrary variation of 77 this 
means that 7; is a solution of the field equation for the 
Lagrangian density L^. □ 



VII. PLANE WAVE SOLUTIONS 

Suppose that M = is Euclidean 3-space equipped 
with Cartesian coordinates x — {x^,x'^,x^) and standard 
Euclidean metric ([B9[) . In this section we construct a 
special class of explicit solutions of the field equations 
for our Lagrangian density ([T9)) . This construction is 
presented in the language of spinors. 

Let us choose Pauli matrices (|B10[) and seek solutions 
of the form 



^{x°,x\x'^,x^) = e-'(P°^''+f=^)C 



(36) 



where po is a real number as in formulas (|5]) and ((D), 
p = (pi,p2,P3) is a real constant covector and 7^ is a 
constant spinor. We shall call solutions of the type ([55)) 
plane wave. In seeking plane wave solutions what we are 
doing is separating out all the variables, namely, the time 
variable x^ and the spatial variables x = {x^,x'^,x^). 

Our dynamic field equation ([?!)) is highly nonlinear so 
it is not a priori clear that one can seek solutions in 
the form of plane waves. However, plane wave solutions 
(1551) are a special case of stationary solutions ([5]) and 
these have already been analyzed in preceding sections. 
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Namely, Theorem [T] gives us an algorithm for the calcu- 
lation of all plane wave solutions by reducing the 
problem to a pair of stationary Weyl equations (O for 
the time- independent spinor field 



(37) 



Substituting formulas (|B2p , (|B10[) and (|37| into equation 
dS]) we get 



TPa +P3 Pi- ip2 
Pi + iP2 TPo - P3 



0. 



(38) 



The determinant of the matrix in the LHS of equation 
([55]) is Pq — p\ — P2~ p\ so this system has a nontriv- 
ial solution C, if and only H p^ — p\ — P2 — pi, — 0. Our 
model is invariant under rotations of the Cartesian coor- 
dinate system (orthogonal transformations of the coor- 
dinate system preserving orientation) so without loss of 
generality we can assume that 



Pi = P2 = 0, P3 = ±po 



(39) 



plane wave solutions we have. To this end, we rewrite 
formula ([^^ in the form 

cos2|po|(a;° + hx^) 
d'^a = ( asin2|po|(a;° + hx^) 

-asin2|po|(a;° + bx^)\ 
d\ = I cos2|po|(a:*' +6x3) 1 ^ = ( I (43) 

where a and h can, independently, take values ±1. It 
may seem that we have a total of 4 different plane wave 
solutions. Recall, however, that we can perform rigid 
rotations of the coframe and that we have agreed (see 
Remark [T] at the end of Section |n| to view coframes that 
differ by a rigid rotation as equivalent. Let us perform a 
rotation of the coordinate system 











hi 


X 


x^ 







where the ± sign is chosen to agree with that in equa- 
tion (j38|) . i.e. upper sign in p9p corresponds to upper 
sign in (I38|) and same for lower signs. Substituting for- 
mulas ([M)) into equation ([55]) and recalling our assump- 
tion dH) we conclude that, up to scaling by a nonzero 
complex factor, we have 



(40) 



simultaneously with a rigid rotation of the coframe 




It is easy to see that the above transformations turn a 
solution of the form P5)) into a solution of this form again 
only with 



Combining formulas (|36p . ((39)) and (|40|) we conclude 
that for each real po 7^ our model admits, up to a rota- 
tion of the coordinate system and complex scaling, two 
plane wave solutions and that these plane wave solutions 
are given by the explicit formula 



a I— > 



-b. 



(41) 



Let us now rewrite the plane wave solutions dH]) in 
terms of our original dynamical variables, coframe § and 
density p. Substituting formulas ((B2|) . (|BT0)) and (|iT|) 
into formulas (|Cl|) - (jC4l) we get p = 1 and 

^ cos2pq{x° ± x^) 
sin2po(a;° ± x^) 


- sin2po(a;'^ ± x^) 

cos2po(x°±a;3) | , z^^^ = | | . (42) 

Note that scaling the spinor C by a nonzero complex fac- 
tor is equivalent to scaling the density p by a positive 
real factor and time shift x^ ^ x^ + const. 

We will now establish how many different (ones that 
cannot be continuously transformed into one another) 



Thus, the numbers a and h on their own do not char- 
acterize different plane wave solutions. Different plane 
wave solutions are characterized by the number c := ab 
which can take two values, -f 1 and —1. 

We have established that for a given positive frequency 
\po I we have two essentially different types of plane wave 
solutions. These can be written, for example, as 

cos2|po|(a;" -I- x^) 
= ( ±sin2|po|(a;°-f 

^Tsin2\po\{x° + x^) 
^\=\ cos2|po|(a;0 + a;3) |, i?3^ = ( | . (44) 

The plane wave solutions (j44p describe traveling waves 
of rotations. Both waves travel with the same velocity 
(speed of light) in the negative x^-direction. The dif- 
ference between the two solutions is in the direction of 
rotation of the coframe: if we fix the spatial coordinate x^ 
and look at the evolution of (|44|) as a function of time x° 
or if we fix time x° and look at the evolution of as 
a function of the spatial coordinate x^ then one solution 
describes a clockwise rotation whereas the other solution 
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describes an anticlockwise rotation. We identify one of 
the solutions ([^^ with a left-handed massless neutrino 
and the other with a right-handed massless antineutrino. 

The bottom line is that our model gives the correct 
number, two, of distinct plane wave solutions. 

VIII. RELATIVISTIC REPRESENTATION OF 
OUR LAGRANGIAN 

In this section we work on the 4-manifold M x M 
equipped with Lorentzian metric ([1]). This manifold is 
an extension of the original 3-manifold M . We use bold 
type for extended quantities. 

We extend our coframe as 

- (ol) ' (45) 



The point of the arguments presented in this section 
was to show that if one adopts the relativistic point of 
view then our Lagrangian density (|19|) takes the espe- 
cially simple form (15(11) . Formula ([50]) is also useful in 
that it allows us to see that our Lagrangian density is in- 
variant under conformal rescalings of the 4-dimensional 
Lorentzian metric g: the arguments from Section |TT] 
(see formulas ([5|)~ PT|) and ((T5)) ) carry over to the 4- 
dimensional setting without change. 

A consistent pursuit of the relativistic approach would 
require the variation of all four elements of the extended 
coframe, giving three extra dynamical degrees of freedom 
(Lorentz boosts in three directions) . We do not do this in 
the current paper, assuming instead that the zeroth ele- 
ment of the extended coframe is specified by formula (j45l) . 

IX. DISCUSSION 



^^«-(^^-J, J = 1,2,3, (46) 

where the bold tensor index ol runs through the values 0, 
1, 2, 3, whereas its non-bold counterpart a runs through 
the values 1, 2, 3. In particular, the Oq in formula P5)) 
stands for a column of three zeros. 

Throughout this section our original 3-dimensional 
coframe § is allowed to depend on time x'^ in an ar- 
bitrary (not necessarily harmonic) manner, as long as 
the kinematic constraint © is maintained. Thus, our 
only restriction on the choice of extended 4-dimensional 
coframe 'd is formula (H5t which says that the zeroth ele- 
ment of the coframe is prescribed as the conormal to the 
original Riemannian 3-manifold M . 

The extended metric ([l} is expressed via the extended 
coframe (^5]) and as 

g = OjkiJj ® -d^ (47) 

where ojk = o^*' := diag(— 1, -|-1) (compare with 

formula (O). The extended axial torsion is 

3 

= ^(-i9"A^rf<+ i9^A -i?^A d^^+'d^A d^^) (48) 

=0 

where d denotes the exterior derivative on R x M (com- 
pare with formula Formula (j48|) can be rewritten as 

rj,ax^yax_^0^^ (49) 

with T'^^ and defined by formulas (|S|) and ([T5|) respec- 
tively. Squaring dH]) we get ||T'^'^||2 = HT'^^lp - ||,9||2 
which implies that our Lagrangian density ^T9\\ can be 
rewritten as 

L(^,p) = ||T-||V. (50) 



The mathematical model presented in Section |TT] is, 
effectively, a special case of the theory of teleparallelism 
Modern reviews of teleparallelism can be found 
m ui. The differences between our mathematical 

model and those commonly used in teleparallelism are as 
follows. 

• We assume the metric to be prescribed (fixed) 
whereas in teleparallelism it is traditional to view 
the metric as a dynamical variable. In other words, 
in teleparallelism it is customary to view (|6]) not 
as a kinematic constraint but as a definition of the 
metric and, consequently, to vary the coframe with- 
out any constraints. This is not surprising as most, 
if not all, authors who contributed to teleparal- 
lelism came to the subject from General Relativity. 

• We take the density of our continuum p to be a 
dynamical variable whereas in teleparallelism the 
tradition is to prescribe it as p = ^/det g . Taking 
p to be a dynamical variable is, of course, equiva- 
lent to introducing an extra real positive scalar field 
p/\/detg into our model 

• We choose a very particular Lagrangian density 
(jSOp containing only one irreducible piece of torsion 
(axial) whereas in teleparallelism it is traditional 
to choose a more general Lagrangian containing all 
three pieces (axial, vector and tensor): see formula 
(26) in In choosing our particular Lagrangian 
density (j50p we were guided by the principles of 
conformal invariance, simplicity and analogy with 
Maxwell's theory. 

The main result of our paper is Theorem[T] which estab- 
lishes that in the stationary setting (prescribed harmonic 
oscillation in time) our mathematical model is equivalent 
to a pair of massless Weyl equations ([2]). The advantage 
of our approach is that it makes the Weyl equation look 
natural to someone with a continuum mechanics back- 
ground. The downside is that our mathematical model 
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is nonlinear which makes it look unnatural to someone 
with a quantum mechanical background. 

The situation here has a certain similarity with in- 
tegrable systems. Say, the Korteweg-de Vries equation 
(mathematical model of waves on shallow water surfaces) 
is nonlinear but the inverse scattering transform reduces 
it to the analysis of a spectral problem for a linear Sturm- 
Liouville operator. In our paper we go the other way 
round, reformulating the spectral problem for the linear 
Weyl operator as a nonlinear equation from continuum 
mechanics. 

From a purely mathematical viewpoint Theorem [1] is 
unusual in that it states that a (particular) second or- 
der partial differential equation is equivalent to a pair of 
first order partial differential equations, which is actually 
hard to believe. Indeed, let us choose a 2-dimensional 
hypersurface S on the 3-manifold M and set a Cauchy 
problem on this surface. When dealing with a second 
order partial differential equation one expects to be able 
to prescribe the value of the spinor field 77 on the surface 
S as well as its normal derivative, whereas when dealing 
with a first order partial differential equation one expects 
to be able to prescribe the value of the spinor field rj only 
(the value of the normal derivative of rj on the surface 
S will be determined by the equation). This argument 
appears to show that there is no way a second order par- 
tial differential equation can be reduced to a pair of first 
order equations. However, our second order partial dif- 
ferential equation happens to be degenerate and does not 
admit the setting of a standard Cauchy problem. This 
degeneracy manifests itself in the property of scaling co- 
variance of our stationary Lagrangian density (|28p . see 
Section IVII for details. Scaling covariance implies that 
our stationary Lagrangian density (I28p vanishes on so- 
lutions of the (second order) field equation which means 
that the value of the spinor field rj on the surface S and 
its normal derivative cannot be chosen independently. In 
order to allay fears that there is something inherently 
wrong with our construction we provide in Appendix |E] 
an elementary example showing by means of an explicit 
calculation that a second order differential equation with 
Lagrangian of the form (I5T]) and ((55t does indeed reduce 
to a pair of first order equations. 

Our construction exhibits a certain similarity with the 
Riccati equation. Recall that the Riccati equation is a 
nonlinear first order differential equation which reduces 
to a linear second order differential equation. We go the 
other way round, reducing a nonlinear second order equa- 
tion to a pair of linear first order equations. However, 
unlike the Riccati equation, our construction works not 
only for ordinary differential equations but also for par- 
tial differential equations. 

Theorem [T] leaves us with two issues unresolved. 

A What can be said about the general case, when the 
spinor field ^ is an arbitrary function of all space- 
time coordinates (x'^ , , x'^ , x^) and is not neces- 
sarily of the form 



B What can be said about the relativistic version of 
our model described in Section IVlIII ? 

The two issues are, of course, related: both arise because 
in formulating our basic model in Section |IT] we adopted 
the Newtonian approach which specifies the time coordi- 
nate x° ( "absolute time" ) . 

We plan to tackle issue A by means of perturbation 
theory. Namely, assuming the metric to be flat (as in 
Section fVIip . we start with a plane wave and then 
seek the unknown spinor field ^ in the form 

^{x°,x\x^,x^) = e~*(P''"''+f•")c(^^a;^x^a;3) (51) 

where C is a slowly varying spinor field. Here "slowly 
varying" means that second derivatives of can be ne- 
glected compared to the first. Our conjecture is that the 
application of a formal perturbation argument will yield 
the Weyl equation ([2]) for the spinor field 

We plan to tackle issue B by means of perturbation 
theory as well. The relativistic version of our model has 
three extra field equations corresponding to the three ex- 
tra dynamical degrees of freedom (Lorentz boosts in three 
directions). Our conjecture is that if we take a solution 
of the nonrelativistic problem which is a perturbation of 
a plane wave (as in the previous paragraph) then, at a 
perturbative level, this solution will automatically satisfy 
the three extra field equations. In other words, we conjec- 
ture that our nonrelativistic model possesses relativistic 
invariance at the perturbative level. 

The detailed analysis of the two issues flagged up above 
will be the subject of a separate paper. 

Appendix A: General notation 

Our general notation mostly follows [l^, the only 
major difference being that we changed the signature of 

Lorentzian metric g^fs from H to — The 

latter is more natural when promoting the Newtonian 
continuum mechanics approach. 

We use Greek letters for tensor (holonomic) indices and 
Latin letters for frame (anholonomic) indices. 

We identify differential forms with covariant antisym- 
metric tensors. Given a pair of real covariant antisym- 
metric tensors P and Q of rank r we define their dot 
product as P • Q ^Po.^...c.^Q h^.p^g"'^' ■ . We 
also define \\Pf := P P. 

All our constructions are local and occur in a neigh- 
borhood of a given point P of the 3-manifold M . We 
allow only changes of local coordinates a;", a — 1,2,3, 
which preserve orientation. 

Working in local coordinates with specified orientation 
allows us to define the Hodge star: we define the action 
of * on a rank r antisymmetric tensor R as 

(*i?)a„+i...o3 := (H)"' v/d^^"'-"'-£ai...a3 (Al) 
where e is the totally antisymmetric quantity, £123 := -|-1. 
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Coframes -d fall into two separate categories, depending 
on the sign of det i^-' q,. We choose to work with coframes 
satisfying the condition 



deti9^c. > 0. 



(A2) 



Condition (jA2p means that orientation encoded in our 
coframe agrees with that encoded in our coordinate sys- 
tem. 

An orthogonal transformation of a coframe is a linear 
map 



-di ^ ^1 = Oh^'" 



(A3) 



where the k are real scalar functions satisfying the con- 
dition SjiO^kO^r = Skr- Of course, orthogonal trans- 
formations map coframes into coframes, i.e. they pre- 
serve the kinematic constraint (jS]). We call an orthog- 
onal transformation special (or a rotation) if the k 
satisfy the additional condition detO-'fc = -1-1. Any two 
coframes satisfying condition (jA2|) are related by a spe- 
cial orthogonal transformation (rotation). 



Appendix B: Spinor notation 

Our spinor notation mostly follows [lB|, the difference 
being that we changed the signature of Lorentzian metric. 

We use two-component complex- valued spinors (Weyl 
spinors) whose indices run through the values 1,2 or 1,2. 
Complex conjugation makes the undotted indices dotted 
and vice versa. 

Define the "metric spinor" 



Cab - - 



-1 

1 



(Bl) 



with the first index enumerating rows and the second 
enumerating columns. We will be using the spinor (|B1[) 
for lowering and raising spinor indices. 
We define 



1 
1 



Qdb 



1 

1 



(B2) 



The spinor (jB2l) can also be used for raising and lowering 
spinor indices. This is a feature of the nonrelativistic 
setting, when we have a specified time coordinate t = x'^ 
and transformations of spatial local coordinates x" , a — 
1, 2, 3, do not involve time. 

Let D be the real vector space of trace-free Hermitian 
2x2 matrices aab ■ Pauli matrices aaab , ol = 1,2,3, are 
a basis in t) satisfying 



(B3) 



where cr^°^ := e^-'^uped^ . Note that formula (jB3|) auto- 
matically implies an analogous formula for the extended 
metric ([1]): 



.... ^ ac I ^ „ ac 
CoLcibf^P +crf3db<^a 



(B4) 



where the bold tensor indices a, (3 run through the values 
0,1,2,3. 

Of course, our Pauli matrices CTq,, a — 1,2,3, are not 
uniquely defined: if ctq, — ffadb are Pauli matrices then so 
are the matrices U*aaU where U is an arbitrary special 
(detU = 1) unitary matrix- function. Note also that un- 
der coordinate transformations our Pauli matrices a : 

aab 

transform as components of a covector: this is indicated 
by the Greek subscript a. 

Let us mention a useful identity for Pauli matrices, 
very similar to (IB4[) but with contraction over tensor in- 
dices instead of spinor ones: 

fctdbCr^cd = — 2eacebd ■ (B5) 
We define the covariant derivatives of spinor fields as 



where P^^j, = r°^b and fi runs through the values 1, 2, 3. 
The explicit formula for the spinor connection coefficients 
r°^f, can be derived from the following two conditions: 



(B6) 



V^a"a6 = 0, (B7) 

where 

ab = C^Cr ab + t p^CT dh — t ^d^ cb ~ ^ ^fcCT ad 

and rl^ay = {a-y} ha'^^'^ag-iS + d^gas " ^a'^q^) are 
the Christoffel symbols. Conditions (jB6p . (IB7p give an 
overdetermined system of linear algebraic equations for 
R-eP^pf), ImP^pfc the unique solution of which is 

1 



-L fib = —^O-a [Oi.L<7 cb~r I fipCF cb) 



(B8) 



Observe that the sign in the RHS of formula (jB8l) is dif- 
ferent from that of formula (A. 9) in [Igj. This is because 
we changed the signature of Lorentzian metric. 
Note that for the standard Euclidean metric 

<?„^ = diag(l,l,l) (B9) 

the traditional choice of Pauli matrices is 

0-ldb = K Qj,cr2d6 = f^ QM,cr3d6 = fQ ^\ (BIO) 



Appendix C: Correspondence between coframes and 
spinors 

In dimension 3 a coframe d and a (positive) density 
p are equivalent to a nonvanishing spinor field ^ modulo 
the sign of ^ in accordance with the formulas 



s = e<^odbe 



(Cl) 
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s^dit^, (C2) 



+ ^ s-^e^Vo^^CVccdC^ (C3) 



= s-leV„af,^^ (C4) 

The above formulas are a special case of those from [13] ■ 
We assume that our Pauli matrices are chosen in such 
a way that the coframe defined by formulas (|C1|) . (|C3|) 
and (|C4|) satisfies condition (jA2p for all 7^ 0- Of course, 
the sign of detz^-'a can always be changed by switching 
from original Pauli matrices to their complex conjugates. 

Note that if we have the standard Euclidean met- 
ric (jB9|) . use traditional Pauli matrices (|B10[) and take 



(C5) 

then formulas ([03]) and (IC4|) give us 

= 5\ . (C6) 



^\ = I -iSe + *<5? 



[curl(t?i+z,92)]„=| 2Vie + ^3ie-^') I- (D3) 
V3(e'+f) , (D4) 



[curl??^]^ = 



doe - do^ 

-2doe 



doe + doe 
-idoe + idoe 





(D5) 



(D6) 



Appendix D: Spinor representation of axial torsion 
and angular velocity 

We show in this appendix that the Hodge dual of axial 
torsion ([8]) is expressed via the spinor field ^ as 



^e^^ocde 



(Dl) 



and that the vector of angular velocity ui defined by for- 
mula (|17l) is expressed via the spinor field ^ as 



^ie^aabdoe - e<yaabdoe) 
e<^ocde 



(D2) 



Note that formulas (jDl[) and (jD2l) are invariant under the 
rescaling of our spinor field by an arbitrary nonvanishing 
real scalar function. 

Formulas (|D1I) and (jD2p are proved by direct substi- 
tution of formulas ([UT]), dSS]) and into © and ^ 
respectively. In order to simplify calculations we observe 
that the expressions in the left- and right-hand sides of 
formulas (|D1I) and (|D2I) have an invariant nature, hence 
it is sufficient to prove these formulas for standard Eu- 
clidean metric (|B9p , traditional Pauli matrices (|B10|) and 
at a point at which the spinor field takes the value (jCSP . 

We have 



,a_ :i+se 
^ - 1 6e 



i+se- se 
i + ise - ise 
^25e 



where curlu :— *du. 

We rewrite the formulas for ^^T^^ and cj in the form 

* j-ax = - id"^) ■ curl(z9i -I- id'^) 

+ 1(^1+ i^2) . curl(i?i - + ■ curl^9^ (D7) 



+ 1(^1 + i^2) X ao(i9^ - ii?^) + ^1?^ X aoz?3 (D8) 

where u ■ v := u^v" (note the absence of complex con- 
jugation) and u X V :— '^{u A v). Substituting formulas 
HI), 03\i and ((D4)) into formula ((D7| we get 

2i 



VaC' + (Vi - iV2)e' - VaC"^ - (Vi + iVa)^' 



which coincides with the RHS of formula (IDip . Substi- 
tuting formulas JCll, ([05]) and (|D6|) into formula ((D8| 
we get 

-*aoe' - idoe 

doe - doe 

which coincides with the RHS of formula (jD2[) . 

An alternative way of proving formulas of the type 
(|Dip and (jD2[) is to choose Pauli matrices CTa, a = 
0, 1, 2, 3, in such a way that a given nonvanishing spinor 
field ^ takes the value (jC5[) in some neighborhood of a 
given point (as opposed to only the point itself). This 
approach was adopted, for example, in [5l. [l8l - [20j . 
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Appendix E: Toy model 

In this appendix we present a toy model showing that 
a second order differential equation with Lagrangian of 
the form (PT|) and (1551) reduces to a pair of first order 
equations. 

We work on the real line R parametrized by the coor- 
dinate X. The dynamical variable (unknown quantity) is 
the scalar function 77 : R — > C \ {0}. Differentiation in x 
is denoted by V. 

Consider a pair of first order linear ordinary differential 
equations 



iVrj ± rj — 0. 
The corresponding Lagrangians are 



(El) 



(E2) 



Equations (jEip are simplified versions of the stationary 
Weyl equations ([5]) and Lagrangians (jE2p are simplified 
versions of the stationary Weyl Lagrangians (|29l) . Note 
that the Lagrangians (jE2[) possess the property of scaling 
covariance ((33|) where : R ^ R is an arbitrary scalar 
function. 

By analogy with ((3T|) . put 



2L+(7?)£_(ry) 
L+irj)-L^i7j) 



i{riVr] — riVr]) 



2^1 



(E3) 



The corresponding field equation (Euler-Lagrange equa- 
tion) is 



./ (Vry) ?7(?7V77-?7V77) rj' 
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2\tj\ 



rj = (E4) 



where the last V in the curly brackets acts on all the 
terms to the right, including those in the square brack- 
ets. Equation (jE4[) is a second order nonlinear ordinary 



differential equation which does not appear to bear any 
resemblance to the first order linear ordinary differential 
equations (|Eip . 

Let us switch to the polar representation of the com- 
plex function rj : 



77 



(E5) 



where r : R ^ (0, +00) and : R ^ R are the new 
dynamical variables (unknown quantities). Substituting 
formula (jE5[) into equation (|E4|) and multiplying by e*"^ 
we arrive at the polar representation of our field equation: 

2i(Vr)(V(^) + r{Vipf + irWip - r = 0. 

Separating the real and imaginary parts we rewrite the 
latter as a system or real equations 



which, in turn, is equivalent to 

V(/3 = Tl, Vr = 0. 



(E6) 



This shows that a complex function 77 is a solution of 
equation (IE4[) if and only if it is a solution of one of the 
two equations (|Eip . 

Of course, the explicit calculations carried out above 
were unnecessary because the toy model considered in 
this appendix is covered by the abstract argument pre- 
sented in Section I VII The point of these explicit calcula- 
tions was to illustrate the degeneracy of field equations 
for Lagrangians of the form ((3T|) and ((33)) : looking at 
(IE6I) one sees the absence of second derivatives. 
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